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SUMMARY 

A modular  model  of  the  sudden-expansion  combustor  has  been 
developed,  combining  a parabolic,  directed-f low  computation 
and  a well-stirred  reactor  model  representing  the  recircula- 
tion region(s).  The  parabolic  module  incorporates  a turbulent 
kinetic  energy  turbulence  model,  and  both  modules  use  the  quasi- 
global  formulation  for  rapid  computation  of  finite-rate  chemical 
reactions.  Coupling  between  the  modules  is  through  a shear 
layer  representation  which  provides  the  parabolic-flow  boundary 
conditions  and  the  stirred-reactor  feed  rates.  Compared  to 
unified  elliptic  models  of  the  dump  combustor,  the  major  advantage 
of  this  approach  is  the  greatly  increased  detail  possible  in 
the  combustor  analysis.  This  detail  in  turn  allows  the  develop- 
ment of  chemical  kinetic  and  flowfield  models  applicable  to 
the  refinement  of  elliptic  formulations.  The  modular  concept 
is  currently  being  applied  in  ramjet  performance  analysis  under 
an  AFAPL  program  and  in  combustor  models  emphasizing  alternate 
fuels  under  DoE  support. 
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ABSTRACT 

This  report  describes  the  development  of  a modular  model 
for  the  prediction  of  the  performance  of  sudden  expansion  burners 
as  a function  of  the  controllable  parameters  relevant  to  comb- 
ustor design.  The  model  is  based  upon  a concept  in  which  the 
recirculation  zone,  treated  as  a stirred  reactor,  is  coupled 
to  a parabolic  boundary  layer  formulation  for  the  flow  outside 
the  recirculation  zone.  Hydrocarbon  oxidation  kinetics  and 
turbulent  kinetic  energy  models  are  employed  in  the  model  develop 
ment.  In  addition  to  the  parabolic-flow  and  stirred  reactor 
elements,  a module  representing  the  fuel  Injection  process  has 
been  developed.  Results  of  the  application  of  the  modular  model 
to  the  analysis  of  cold-flow  and  reacting-flow  dump  combustor 
experimental  data  are  described. 
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NOMENCLATURE  (coat.) 


MEANING  EQUATION 

temperature  12 , 18 

axial  mean  velocity  component  13 

axial  mean  velocity  component 
on  centerline  23 

axial  mean  velocity  component  on 
centerline  at  X/D  =0  23 

axial  turbulent  velocity  component  20 

radial  mean  velocity  component  13 

stirred  reactor  volume  7 

radial  turbulent  velocity  component  20 

azimuthal  turbulent  velocity  component  20 
chemical  production  term  for  species 
i 11,14 

production  of  kth  gas  specie  due 
to  homogeneous  gas  phase  reactions  7 

production  of  kth  particle  class  7 

axial  coordinate  13 

lateral  coordinate  in  parabolic 
mixing  model  13 

mass  fraction  of  species  1 7,14 

eddy  kinematic  viscosity  10,11 

eddy  kinematic  viscosity  at  half- 
velocity point  25 

gas  phase  eddy  viscosity  14 

gas  thermal  conductivity  10 

turbulent  dynamic  viscosity,  = pe^  19 


April  1978 


SAI-78-008-WH 


NOMENCLATURE  CCout.) 

SYMBOL 

MEANING 

EQUATION 

P 

density 

12 

characteristic  density  of 

recirculation  region 

10,11 

Tt 

shear  stress 

22 

stream  function 

- 

SUBSCRIPTS 

(>i 

species  i 

P 

outside  dividing  streamline 

R 

inside  dividing  streamline 

w 

dividing  streamline 

SUPERSCRIPTS 

o' 

Influx  to  stirred  reactor 

outflow  from  stirred  reactor 


April  1978 


SAI-78-008-WH 


I . INTRODUCTION 

Current  and  future  airbreathing  propulsion  system  concepts 
are  tending  toward  lower  weight  and  smaller  volume  designs. 

An  example  is  the  Integral  rocket/ramjet  engine  which  uses  a 
common  combustor  for  rocket  boost  and  sustained  ramjet  flight. 
Unfortunately,  the  current  design  trends  create  stringent  re- 
quirements for  higher  performance  from  smaller  volumes,  and 
the  combustion  processes  are  not  well  enough  understood  so  that 
combustor  development  has  been  carried  out  largely  by  employing 
costly  cut-and-try  methods. 

The  basic  difficulty  in  understanding  combustor  behavior 
is  the  coupling  between  turbulent  mixing  and  reaction  kinetics. 
This  problem  is  not  unique  to  combustor  flows,  but  is  encountered 
in  other  systems  where  flow,  mixing,  and  reaction  times  are 
all  of  the  same  order  of  magnitude,  for  example,  external  burn- 
ing devices  for  thrust  and/or  control  and  in  gas  dynamic  and 
chemical  laser  systems.  In  addition,  many  of  the  reacting 
flows  of  general  interest  Involve  condensed  phase  species  and 
large  scale  recirculation  zones  such  as  are  found  in  the  "dump" 
burner  concept.  Problems  in  ignition,  flame  stabilization, 
combustion  efficiency,  and  flowfield  and  combustion  instabilities 
also  occur  in  the  dump  burner  as  well  as  in  many  other  combustion 
systems  of  Interest,  and  these  problems  need  to  be  better  under- 
stood to  support  engineering  Judgement. 

Consequently,  continued  research  is  required  in  both  sub- 
sonic and  supersonic  flows  in  order  to  gain  the  necessary  under- 
standing of  the  controlling  mechanisms.  Only  in  this  way  can 
effective  design  principles  be  developed  to  provide  optimum 
propulsion  system  performance  over  wide  operating  ranges. 

The  objective  of  the  work  summarized  in  this  report  is  to 
eliminate  the  deficiencies  Just  described  and  provide  a sound 
basis  for  effective  combustor  design  criteria,  through  the 
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development  of  analytical  models  for  the  sudden-expansion  comb- 
ustor flowfield.  The  major  emphasis  in  this  program  is  on  the 
development  of  a "modular"  model  for  the  sudden  expansion  comb- 
ustor, in  which  proved  computatioral  elements,  including  a 
parabolic  flowfield  analysis  procedure  and  a perfectly-stirred 
reactor  formulation  are  coupled  together  to  model  the  overall 
flowfield.  In  this  report  the  features  of  the  modular  approach 
are  reported  and  progress  in  the  development  and  application 
of  the  model  to  cold-flow  and  reacting-flow  combustor  data  is 
described.  This  work  has  involved,  and  will  continue  to  in- 
volve, close  coordination  with  the  experimental  work  in  progress 
at  the  Air  Force  Aero  Propulsion  Laboratories  at  Wright-Patterson 
Air  Force  Base,  as  well  as  other  related  programs. 
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II.  MODULAR  MODEL  CONCEPT 

The  basic  characteristics  of  the  sudden  expansion  burner 
are  depicted  in  Figure  1.  In  essence,  the  key  feature  of  the 
flow  field  is  the  recirculation  region  just  downstream  of  the 
expansion,  separated  from  the  directed  flow  by  the  dividing 
streamline.  Because  of  this  recirculation  region,  in  which 
there  is  no  characteristic  flow  direction,  the  flow  field  in 
the  sudden  expansion  burner  is  in  general  described  by  ellip- 
tic equations;  i.e.,  the  boundary  layer  approximations  cannot 
be  applied  to  the  whole  flow  and  the  elliptic  Navier-Stokes 
equations  must  be  used.  However,  it  is  possible,  at  least 
conceptually,  to  divide  the  flow  field  into  two  regions.  The 
first,  the  directed  flow  region,  consists  of  that  part  of  the 
flow  between  the  dividing  streamline  and  the  burner  centerline. 

The  second  region  is  the  recirculating  flow  region,  and  the  two 
regions  are  coupled  through  the  conditions  at  the  dividing 
streamline. 

The  concept  of  the  modular  approach  to  the  computation  of 
the  flow  in  a sudden  expansion  burner  involves  the  application 
of  available  techniques  to  the  separate  regions  of  the  flow 
field.  Thus,  the  directed  flow  region  is  treated  using  a para- 
bolic finite  difference  solution  of  the  boundary  layer  form  of 
the  equations  of  motion.  The  recirculation  regions  are  treated 
as  perfectly  stirred  reactors,  with  the  reactor  feed  rate 
defined  by  the  species  diffusion  flux  across  the  dividing  stream- 
line. The  shape  of  the  dividing  streamline  must  be  specified 
a priori;  the  resulting  wall  static  pressure  distribution  may 
then  be  compared  with  experiment  and  another  dividing  streamline 
shape  specified,  if  necessary.  Since  the  flow  field  is  in  general 
elliptic,  the  modular  approach  approximates  the  true  flow  field, 
and  the  specification  of  the  species  fluxes  across  the  dividing 
streamline  must  be  carried  out  iteratively.  However,  despite 
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Figure  1.  Schematic  of  Sudden  Expansion  (Dump) 
Burner 
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its  limitations,  this  method  does  allow  the  application  of 
well-tested  and  proved  flow  models  to  the  sudden  expansion 
burner  geometry,  and  is  capable  of  yielding  parametric  infor- 
mation of  Interest  in  sudden  expansion  burner  design. 
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III.  MODULAR  MODEL  FORMULATION 

The  basic  elements  of  the  current  modular  model  formulation 
are  the  parabolic  finite-difference  computational  technique 
developed  by  Boccio,  Weilerstein  and  Edelman  [ l] , modified  to 
incorporate  the  turbulent  kinetic  energy  model  developed  by 
Harsha  [2],  used  for  the  directed-f low  portion  of  the  analysis, 
and  the  stirred  reactor  computation  developed  by  Edelman  and 
Weilerstein  [ 3] . Both  of  these  elements  make  use  of  the  quasl- 
global  model  developed  by  Edelman  and  Fortune  w for  rapid 
computation  of  finite-rate  hydrocarbon  air  kinetics.  These 
elements,  or  modules,  are  coupled  together  through  a simplified 
representation  of  the  turbulent  shear  layer  which  exists  between 
the  directed  flow  and  the  recirculation  region.  In  the  modular 
approach,  the  shear  layer  representation  is  used  to  define  the 
gradients  in  velocity,  species,  and  enthalpy  between  the  two 
regions  of  the  flow,  thus  providing  both  the  boundary  conditions 
on  the  directed  flow  and  the  stirred  reactor  feed  rates.  Details 
of  the  mathematical  formulations  used  for  the  flowfield  region 
are  discussed  in  this  section;  details  of  the  numerical  analysis 
procedure  may  be  found  in  Ref . [ 5] . 

1.  Coupling  Relations:  The  Shear  Layer  Model 

The  key  feature  of  the  modular  concept  lies  in  the  coupling 
relations  asswned  along  the  dividing  streamline  that  (1)  supply 
the  boundary  condition  for  the  parabolic  computations  and  (2) 
determine  the  feed  rates  for  the  perfectly  stirred  reactor 
computation.  These  coupling  conditions  are  obtained  through 
a simplified  model  of  the  turbulent  shear  layer  separating  the 
recirculation  zone  from  the  directed  flow,  as  is  shown  schemat- 
ically in  Figure  2.  Note  in  the  following  that  the  shear  layer 
is  assumed  to  be  vanishingly  thin,  that  is,  the  gradients 
established  through  the  shear  layer  model  are  applied  along  the 
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dividing  streamline. 


X 


Figure  2.  Coupling  Relation  Definition  Sketch 

It  is  assumed  that  the  shear  layer  can  be  modeled  as  a region 
of  width  £(x)  across  which  all  dependent  variables  (i.e.,  velocity, 
temperature,  turbulent  kinetic  energy,  and  species  mass  fractions) 
vary  linearly.  The  width  of  the  shear  layer  is  itself  assumed 
to  be  specifiable  by  the  relation 

i=a+bx  (1) 

in  which  a and  b are  constants.  Since,  for  planar  shear  layers, 
the  growth  rate  can  be  related  to  the  dividing  streamline  shear 
stress,  specification  of  the  constant  b Implies  a shear  stress 
distribution  along  the  dividing  streamline,  while  the  constant 
a can  be  taken  to  be  a parameter  which  is  related  to  the  boundary 
layer  thickness  at  the  step. 

As  shown  in  Figure  2,  the  dividing  streamline  shape,  R„(x) 
is  specified;  this  specification  also  defines  the  dividing 
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streamline  stream  function  value  Along  the  turbulent 

shear  stress  is  defined  as  described  above,  and  thus  the 
turbulent  kinetic  energy  is  defined  through  the  relation 

= 0.3  pk^  (2) 

To  obtain  the  gradients  in  temperature  and  species  mass  fraction, 
linear  profiles  across  the  shear  layer  are  assumed.  If  is 
defined  as  the  streamline  immediately  "outside"  the  recirculation 
region  (i.e.,  within  the  directed  flow)  and  is  the  streamline 
immediately  "inside"  the  recirculation  zone,  then,  for  species 
i 


and  for  the  temperature 


where  and  Tj^  are  the  values  of  species  mass  fractions  and 
temperature  obtained  from  the  stirred  reactor  solution. 

The  overall  flowfield  computation  using  the  modular  approach 
proceeds  as  follows.  A dividing  streamline  shape  R (x)  is  assumed, 
and  the  shear  layer  width  expression  and  shear  stress  distribution 
is  defined.  An  initial  state  for  the  stirred  reactor  computation 
is  assumed,  resulting  in  values  of  and  Tj^.  With  these  data 
available  as  boundary  conditions,  the  parabolic  mixing  calculation 
is  carried  out  to  the  end  of  the  recirculation  zone.  This  calcula- 
tion defines  the  species  mass  fraction  and  temperature  gradients 
at  iLDd  these  values  are  used  to  obtain  the  stirred  reactor 

feed  rates.  A new  stirred  reactor  computation  is  carried  out 
using  the  new  feed  rates,  resulting  in  a new  specification  of 
and  Tj^,  and  the  parabolic  computation  is  repeated.  The 
procedure  is  repeated  until  the  changes  in  the  stirred  reactor 
composition  from  computation  to  computation  become  small,  at 
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which  point  the  coupling  iteration  has  converged.  The  parabolic 
calculation  is  then  carried  out  to  the  end  of  the  combustion 
chamber,  completing  the  solution. 

2.  Well-Stirred  Reactor:  The  Recirculation  Zone  Model 

Flowfield  regions  in  which  intense  backmixing  occurs  can 
approach  the  limit  of  complete  mixing,  and  thus  the  well-stirred 
reactor  concept  is  attractive  for  representing  the  recirculation 
region  in  the  sudden  expansion  burner.  Because  the  mixing  process 
is  considered  to  be  complete,  the  mathematical  formulation  for 
the  perfectly  stirred  reactor  is  straightforward.  The  describing 
equations  are  given  as  follows: 


Global  Conservation  of  Mass; 


rb  = constant  - Z *?  = Z (5) 

k k k 

where  k is  the  k^^  component  which  can  include  i gas-phase 
species  and  j droplet  types.  A droplet  type  (or  class)  is 
defined  by  its  composition  (fuel,  water,  etc.)  and  its  size. 

The  superscripts  I and  O refer  to  inflow  and  outflow,  respectively. 


Conservation  of  Energy: 
h = h^  + Q/m 

where  Q is  the  net  rate  of  external  heat  addition  to  the  zone 
and  the  inflow  of  enthalpy  is  written  to  allow  each  component 
to  enter  the  reactor  with  an  arbitrary  temperature. 


(6) 


i 


Conservation  of  Species: 

The  conservation  of  species  k requires  that  its  net  outflow 
be  equal  to  the  rate  of  production  due  to  chemical  and  phase 
transition  processes.  Formally,  this  principle  gives  k 
equations  of  the  form; 


‘k 


“k  * 


(7) 
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where  superscripts  G and  P refer  to  homogeneous  gas-phase  and 
dropwlse  production  rates,  respectively. 

Although  there  exist  more  or  less  standard  techniques  for 
the  solution  of  such  algebraic  equations,  an  alternative  ap- 
proach Involves  the  use  of  the  non-steady  form  of  the  species 
conservation  equations,  viz: 


da 


k 


dt 


m 


/w  /i  \3 

A - 


(8) 


where  t is  the  time  variable  of  significance  only  during  the 
transient  period.  Thus,  Equation  (8)  is  identical  to  Equation 
(7)  when: 


da 

dt 


K 


(9) 


The  remaining  working  equations  retain  their  steady  state  form. 

These  equations,  along  with  expressions  for  the  volumetric 
production  rates  appearing  in  equation  (8) ,expresslons  for  the 
enthalpy  as  a function  of  species  concentrations  and  the  equation 
of  state,  define  the  temperature  and  species  concentrations 
in  the  stirred  reactor,  given  the  inflow  rates  for  species  and 
enthalpy.  In  the  modular  model,  the  net  inflow  of  species  and 
enthalpy  can  be  expressed  as  line  integrals  involving  gradients 
evaluated  along  the  dividing  streamline,  so  that  for  the  modular 
model  the  energy  and  species  conservation  equations  for  the 
stirred  reactor  can  be  written 
ENERGY 


-2 


ds 


/■®  V 

''To  3®?  1 

- 2"  y Rc<*>p‘d  4 

L^i  7F"  J 

(10) 

10 


April  1978 


SAI-78-008-WH 


SPECIES 


+ 


neglecting  the  dropwise  species  production  rate  term.  Here 
the  superscript  I refers  to  inflow  into  the  stirred  reactor 
(recirculation  region)  and  superscript  o refers  to  outflow 
from  the  stirred  reactor,  V is  the  reactor  volume,  p a character- 
istic  density  of  the  stirred  reactor  region,  and  pc^  a character- 
istic eddy  diffusivity,  evaluated  from  the  outer  flow  field 
solution  in  the  region  of  the  dividing  streamline.  The  term 
represents  the  rate  of  production  of  species  i caused  by 
chemical  reactions  and  Q represents  the  heat  input  to  the 
stirred  reactor  region  through  the  combustor  walls. 

Equations  (10)  and  (11),  along  with  an  equation  of  state 


P 


where 

R - the  universal  gas  constant 

= the  molecular  weight  of  species  i 
are  used  to  establish  a new  reactor  state,  given  the  feed  rates 
established  from  the  parabolic  directed  flow  solution;  i.e., 
given  the  diffusive  fluxes  of  species  and  energy  across  the 
dividing  streamline.  Equations  (10),  (11)  and  (12)  are  solved 
for  new  values  of  the  and  the  recirculation  zone  temperature, 
Tj^.  The  process  is  then  iterated  until  changes  in  the  para- 
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meters  describing  the  state  of  the  recirculation  zone  are  small. 

3.  Parabolic  Mixing:  The  Directed-Flow  Model 

The  next  major  element  of  the  modular  model  for  a sudden 
expansion  combustor  is  the  formulation  for  the  directed  flow 
portion  of  the  combustor  flowfield.  It  is  assumed  that  the 
boundary  layer  approximations  apply  to  this  part  of  the  flow- 
field,  so  that  the  describing  equations  are  parabolic.  For  a 
steady,  axisymmetrlc  flow  these  equations  may  be  written: 

Global  Continutly 


3ypu  + 3ypv  ^ 0 
dx  3y 


Species  Diffustion  for  the  i Specie: 

^“i  . ®“i  1 3 j 

ai-  * P''  ^ ? W i ' 

Uomentum  Equation: 


3u  . 3u  1 1 3 X 3u  X 

35E  ^ " 7 j ly  ^ ^ 


(13) 


( £ 

rsa.-l 

1 

y P ^ 

Lsy 

( c 

(14) 


^ x)  ^ 

3y  3x 


(15) 


and  the  energy  equation 


1 ^ * <V1)  ly(#)|  ' 

in  which  P^.  and  S^  represent  the  Prandtl  and  Schmidt  numbers, 
respectively.  These  equations,  along  with  expressions  for  the 


(16) 
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enthalpy 


H « h + u^/2  I 
h = I h^(T)  j 

and  the  equation  of  state 
p = pRT  E (o^/M^) 


(17) 


(18) 


can  be  solved,  given  an  expression  for  the  turbulent  momentum 
diffusivity  e^,  or  the  turbulent  eddy  viscosity  y,p  = 

The  problem  of  the  most  appropriate  formulation  for  the  eddy 
viscosity  in  a turbulent  flow  has  occupied  the  attention  of 
researchers  in  turbulent  flow  for  many  years,  and  numerous 
proposals  for  the  appropriate  form  of  the  eddy  viscosity  have 
been  made.  In  general,  up  until  a few  years  ago,  models  for 
the  turbulent  eddy  viscosity  involved  a relatji-on  between  a 
local  length  scale  and  a local  measure  of  the  velocity  gradient. 
The  free  turbulent  mixing  model  proposed  by  Prandtl  in  1942  [6] 
is  a case  in  point;  with  this  model 

y,p  = pcb|Au|  (19) 

where  c is  a constant,  b is  a measure  of  the  width  of  the  mix- 
ing region,  and  Au  a measure  of  the  velocity  difference  across 
the  mixing  region. 

While  eddy  viscosity  models  have  in  certain  circumstances 
enabled  successful  calculations  of  particular  turbulent  flows 
to  be  carried  out,  in  general,  the  constant  of  proportionality 
involved  in  the  model  varies  markedly  and  unpredictably  in 
different  flow  fields  f7] . 

The  observed  lack  of  generality  of  eddy  viscosity  models 
has  led  to  the  development  of  a class  of  models  in  which  addi- 
tional partial  differential  equations  are  written  to  obtain 
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the  spatial  variation  of  the  turbulent  shear  stress.  Of  these, 
the  most  highly  developed  are  the  turbulent  kinetic  energy 
models  (8,  9). 

If  the  turbulent  kinetic  energy  is  defined  as 

k « 1/2  (u’^  + v'^  + w’^)  (20) 


an  equation  describing  the  spatial  variation  of  this  quantity 
can  be  obtained  from  the  momentum  and  continuity  equations  by 
appropriate  manipulations. 

For  steady  axisymmetric  compressible  flow  the  turbulent 
kinetic  energy  equation  can  be  written  [2] 


pu  3k  pv3k  _ 1 _3  / 

3x  3y  ~ y 3y 


,luf 


IV 


(21) 


in  which  the  various  terms  are  interpreted  as: 

I transport  by  mean  flow  (advection) 

II  diffusion 

III  production 

IV  dissipation 

In  order  to  use  Equation  (21),  a relationship  between  the  turb- 
ulent shear  stress,  x,  and  the  turbulent  kinetic  energy,  k, 
must  be  postulated.  One  such  relationship  is 

Tt  = aj^pk  (22) 

For  most  turbulent  shear  flows  the  parameter  a^  must  be  allowed 
to  vary  laterally.  In  general,  for  an  axisymmetric  jet  or 
ducted  flow,  a^  can  be  written  [2] 
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for  Og  >0.9 


for  U < 0.9  0.  and  y < y max 
c — j 


(23) 


(24) 


for  0 < 0.9  O4  and  y > y max 
c “ j 


where  U represents  the  local  axial  velocity  component  on  the 

C 

centerline,  Uj  the  value  of  the  centerline  velocity  at  the  jet 


origin,  y__„  the  value  of  y at  which 

lUttX 


3u 

3y 


3u 

3y 


the  maximum  value  of 


max 


3u 

W • 


is  a maximum,  and 
Further  details  regarding 


the  lateral  variation  of  a^^  may  be  found  in  Reference  2. 

To  describe  the  variation  of  ag,  which  appears  in  the 
dissipation  term  of  Equation  (21),  a "turbulent  Reynolds  number" 
defined  by 


Aul, 


m 


(25) 


This  parameter  is  a 


is  used  as  the  independent  variable, 
function  of  x only, 
viscosity  at  the  maximum  shear  point  in  a lateral  profile,  i.e 


The  term  represents  the  eddy  kinematic 


m 


(26) 


and  is  a characteristic  length  scale,  effectively  the  lateral 
width  of  the  turbulent  flow.  This  scale,  which  is  also  used 
in  the  dissipation  expression  Equation  (21),  is  determined  by  the 
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following  set  of  relationships. 

For  the  two-dimensional  shear  layer,  or  first  regime  of 
a jet  (i.e.,  for  = Uj) 

where  Ay  is  the  lateral  distance  between  the  points  at  which 
(u-u  )/(u.-u  ) = 0.95  and  (u-u  )/(u.-u  ) = 0.05,  for  profiles 


(27) 


-e^/ x-j  -e'  e- 

which  are  not  fully  developed.  In  these  expressions,  u 


IS 


the  velocity  at  the  outer  edge  of  a jet;  at  a wall  u^  = 0. 
For  fully  developed  profiles,  which  can  be  approximated  by  a 
cosine  function. 


= 1.57  Au/(3u/3y) 


m 


(28) 


where  the  subscript  m again  refers  to  the  position  in  a lateral 
profile  at  which  the  shear  stress  is  a maximum.  Equation  (28) 
provides  a characteristic  length  scale  determination  which 
avoids  the  problem  of  defining  the  edge  of  an  asymptotic  pro- 
file. Because  the  cosine  function  is  not  always  a good  approx- 
imation to  the  true  computed  velocity  profile,  in  practice 
Equation  (28)  is  only  used  in  place  of  Equation  (27)  wherever 
the  value  of  determined  by  it  satisfies  the  inequality 
Ay  < ili,<  1.57  Ay. 


1 "k- 

For  the  second  regime  of  jets,  i.e. 


^c  ^ ^j 


when 


*'k  “ ^^1/2 

where  value  of  r at  which  (u-Ug)/(u^-Ug)  = 0.5, 


(29) 


In  the  definition  of  R„ 


Au  = u.-u  in  the  first  regime 
J ® 


of  jets,  and  Au  = u-u  in  the  second  regime  of  jets. 

The  value  of  the  dissipation  constant,  ag,  is  related  to 
the  local  value  of  R,j,  and  to  the  initial  jet  density  ratio. 
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These  relationships  were  obtained  through  comparison  of  pre- 
diction with  a wide  variety  of  jet  and  wake  experiments  [2]; 
in  summary,  they  are: 

Variation  in  ^2  ''i** 


for  0 < R,p  < 185, 
for  185  < < 360/ 

for  > 360, 


^2  “ 1-65 

a.2  ~ 0.46  + 0.00762R^ 
a2  - 3.20 


Variation  in  a2  with  density  ratio 


For  nonunity  initial  density  ratio,  multiply  a2  by  1/c 
where , for 


Pel 

> 1 c,  = 0.984  + 0.016  - — 
Pji  1 “ji 


and  for 


^ < 1 - 0.9S  + O.OS^R^T^^j  / (Vtj) 


April  1978 


SAI-78-008-WH 


In  these  relations,  the  subscript  1 refers  to  conditions  at 
the  origin  of  mixing,  R is  the  gas  constant  appropriate  to 
the  outer  flow,  and  T,p^  is  the  total  temperature  of  that  flow. 
Equation  (32)  is  used  when  the  Jet  density  is  greater  than  the 
free-stream  density  in  order  to  preserve  the  ^2'^  relationship 
for  the  supersonic  shear-layer  case,  for  which  the  "Jet"  side 
is  taken  to  be  the  high-speed  stream.  Thus,  for  the  case  of 
compressible  but  identical-gas  mixing  problems,  the  value  of 
c^  obtained  through  use  of  Equation  (32)  is  unity,  while  for 
gases  with  different  molecular  weights  or  stagnation  temperatures 
a nonunity  c^  is  obtained. 

The  turbulent  kinetic  energy  model  described  in  this  section 
has  been  successfully  applied  to  a wide  variety  of  flows,  with 
and  without  variation  in  density.  For  free  flows,  numerous 
examples  of  the  applications  of  this  model  are  given  in  Ref- 
erence 2;  it  has  also  been  applied  recently  to  the  computation 
of  ducted  flows  [lO]. 

The  complete  set  of  describing  equations  for  the  parabolic 
flowfield  is  thus  made  up  of  equations  (13)  - (18)  along  with 
equation  (21)  and  its  associated  algebraic  relations.  With 
boundary  conditions  prescribed  from  the  stirred  reactor  solution 
through  use  of  the  coupling  conditions  described  in  Section 
III.l,  these  equations  are  numerically  integrated  using  the 
solution  technique  of  Reference  [1]  as  described  in  detail 
in  Reference  [s]. 

4.  Chemical  Kinetics:  The  Quasi-Global  Model 

In  both  the  stirred  reactor  and  directed  flow  portions  of 
the  modular  model  the  volumetric  production  rate  terms  appearing 
in  the  species  transport  equations  are  evaluated  using  a full 
hydrocarbon  chemical  kinetics  scheme  based  on  the  quasi-global 
kinetics  model  [4].  This  model  has  as  a key  element  a subglobal 
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oxidation  step 

* f O2  i ^2  * <33) 

This  reaction  is  unidirectional  with  an  empirically  determined 
rate  (grams  of  fuel/cc/sec)  given  by 

A tV-3  [^n®m]*  [^2]  ®*P[-  (E/BT)]  (34) 

with  the  constants  A,  b,  and  E/R  defined  in  Table  1,  where  P 
must  be  given  in  atmospheres,  T in  degrees  Kelvin  and  [ ] 
denotes  molar  concentration. 

Coupled  to  this  subglobal  step  are  the  intermediate 
reversible  reactions  given  in  Table  1.  Some  observations  on 
this  scheme  should  be  noted.  First,  it  is  generally  agreed 
that  in  hydrogen  bearing  systems  carbon  monoxide  is  most  rapidly 
oxidized  to  CO2  via  the  reaction 

CO  + OH  ^ COg  + H (35) 

We  have,  however,  included  other  reactions  involving  CO  and 
oxygen  for  completness  Including 

CO  + Og  ^ COg  +0  (36) 

and 

CO  + O + M ^ COg  + H (37) 

where  M is  the  general  third  body.  Reactions  (36)  and  (37)  are 
much  slower  than  reaction  (35),  but  their  inclusion  was  necessary 
for  basic  studies  performed  on  systems  including  the  CO/air  system. 
In  addition,  a number  of  reactions  involving  NO  are  included  which 
represent  a necessary  extentsion  of  the  basic  Zeldovich  mechanism 
to  account  for  certain  of  the  ambient  long  time  NO-to-NOg  con- 
version reactions  which  occur  in  the  atmosphere  particularly 
when  coupled  with  appropriate  daylight  photochemical  mechanisms. 
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Table  1.  lataaM  C-l>0  ehaalcal  ktaatle  rcaetloa  aacbaaln  aapC-^/IT) 
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IV.  MODELS  FOR  THE  FUEL  INJECTION  PROCESS 

In  addition  to  the  basic  components  of  a parabolic,  directed 
flow  analysis  and  a well-stirred  reactor  formulation,  the  modular 
concept  can  be  extended  to  include  modules  which  represent  other 
elements  of  the  dump  combustor  flowfield,  for  example,  the  fuel 
injection  process.  The  detail  of  the  computation  provided  by 
the  use  of  a parabolic  directed  flow  analysis  (as  opposed  to 
the  relative  coarseness  of  the  nvimerical  grid  allowable  in 
current,  unified,  elliptic  solution  techniques)  is  the  key  feature 
of  the  modular  model  that  allows  the  Inclusion  of  a fuel  in- 
jection module  in  the  complete  analysis.  This  is  particularly 
true  in  the  case  of  liquid  fuel  injection,  for  at  the  fuel/air 
ratios  appropriate  for  dump  combustor  operation,  the  liquid 
fuel  streams  initially  occupy  a very  small  portion  of  the  over- 
all combustor  cross-sectional  area.  It  is  also  worth  noting 
that  it  has  been  stated  that  the  further  development  of  unified, 
elliptic  combustor  flowfield  calculations  requires  the  develop- 
ment of  what  are  essentially  modular  models  for  such  features 
as  the  fuel  injection  process  (Reference  11)  since  these  key 
processes  occur  on  scales  much  smaller  that  feasible  numerical 
resolution  allows.  Thus  the  development  of  fuel  injection, 
vaporization,  and  spreading  models  through  use  of  the  overall 
modular  model  concept  is  necessary  for  further  development  of 
unified  techniques. 

The  liquid  fuel  injection  model  makes  use  of  a combination 
of  empirical  information  and  turbulent  mixing  calculations. 

For  example,  the  fuel  jet  penetration  from  the  wall  is  computed 
through  the  use  of  the  penetration  correlation  developed  by 
Catton,  Hill  and  McRae  (Ref.  12),  using  the  breakup  time  correla- 
tion developed  by  Clark  (Ref.  13)  to  compute  the  downstream 
distance  at  which  penetration  is  to  be  computed.  That  is,  it 
is  assvuned  that  the  fuel  jet  has  turned  paralled  to  the  airflow 
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at  the  axial  position  at  which  the  inital  fuel  Jet  has  broken 
[ up  into  droplets,  as  given  by  the  breakup  time  correlation  and 

I the  local  airflow  velocity. 

Since  the  basic  modular  model  formulation  involves  an 
axisymmetric  flowfield,  individual  fuel  jets  cannot  be  resolved, 
and  it  is  assumed  in  the  formulation  that  the  liquid  fuel  spray 
forms  an  annulus  whose  cross-sectional  area  may  either  be 
specified  or  computed  based  on  an  assumed  fuel  spray  bulk 
velocity.  A bulk  spray  evaporation  correlation  is  then  used 
to  compute  the  fuel  vaporization  rate;  this  correlation,  developed 
by  Ingebo  and  Foster  [14]  is  a function  of  the  initial  velocity 
and  temperature  difference  between  the  fuel  spray  and  the 
surrounding  air  stream.  Spreading  of  the  fuel  jet  is  computed 
through  use  of  a turbulent  mixing  hypothesis  as  for  the  mixing 
■ process  in  the  remainder  of  the  parabolic  flow. 

Figure  3 shows  the  results  of  a computation  of  the  fuel 
' injection  process  for  three  fuel  injectors,  located  in  the 

combustor  inlet  wall,  along  the  centerline,  and  in  a midstream 
position.  Shown  are  the  computed  contours  of  the  fuel  mass 
i fraction,  Up,  with  the  vapor-phase  fuel  shown  as  the  solid 

line  and  the  liquid  phase  fuel  as  the  dotted  line,  as  a function 
of  both  axial  and  radial  position  in  the  combustor  inlet.  For 
these  calculations  the  air  inlet  velocity  was  approximately 
700  ft/sec  at  a temperature  of  1600  °K;  the  overall  fuel-air 
equivalence  ratio  was  0.6.  The  results  shown  In  Figure  3 provide 
a good  example  of  the  detail  of  the  fuel  Injection  process 
available  through  use  of  this  aspect  of  modular  modeling. 
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AXIAL  POSITION,  FT 

Figure  3.  Predicted  Fuel  Mass  Fraction  Distributions  for 

Simultaneous  Injection  at  Three  Radial  Location. 
Inlet  Air  Velocity  700  ft/sec.  Inlet  Air 
Temperature  1600  ^K,  Overall  Equivalence  Ratio  0.6 
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V.  RESULTS  OF  APPLICATION  OF  THE  MODULAR  MODEL 


During  this  contract  period,  development  work  on  the  modular 
diunp  combustor  analysis  has  centered  on  combining  the  elements 
of  the  analysis  and  determining  the  sensitivity  of  the  results 
obtained  using  the  approach  to  the  assumptions  involved  in  the 
model.  The  basic  accuracy  of  the  directed  flow  module,  in- 
corporating a finite-rate  kinetics  mechanism  and  the  one- 
equation  TKE  turbulence  model,  was  earlier  tested  in  a series 
of  computations  of  a reacting  hydrogen-air  Jet,  as  reported  in 
Ref.  15. 

In  order  to  test  the  predictions  of  the  modular  approach, 
a series  of  calculations  were  carried  out  for  comparison  with 
the  cold-flow  experimental  results  obtained  by  Drewry  (Ref.  16). 

The  configuration  tested  in  these  experiments  is  shown  schematically 
in  Figure  4,  which  also  gives  the  initial  conditions  for  the 
experiment  and  computations.  Note  that  the  length  of  the 
recirculation  zone  used  in  the  calculation  is  that  measured 
by  Drewry,  as  shown  in  Figure  4. 

The  sensitivity  of  the  computed  wall  static  pressure  pro- 
file to  the  assumptions  made  in  the  model  about  the  shape  of 
the  dividing  streamline  and  the  dividing  streamline  shear  stress 
distribution  is  shown  in  Figure  4.  Initially,  the  dividing 
streamline  was  assumed  to  have  a linear  shape,  and  a shear  stress 
level  typical  of  a pipe  flow  was  chosen.  The  shear  stress  is 
represented  on  Figure  4 as  a skin  friction  coefficient,  where 

for  ease  in  comparison  of  the  shear  stress  levels  with  values 
appropriate  for  ducted  flows.  It  can  be  seen  that  the  pipe- 
flow  value  of  produces  a considerably  larger  static  pressure 
rise  than  is  experimentally  observed,  and  that  the  assumption 
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of  a linear  shape  produces  a too  rapid  rate  of  increase  of 
pressure.  Thus,  the  shape  of  the  dividing  streamline  was 
changed  to  a parabola,  such  that  the  overall  length  of  the 
recirculation  region  was  unchanged  and  the  dividing  streamline 
was  initially  parallel  to  the  ccunbustor  wall  at  the  dump  station. 
As  can  be  seen  from  Fig.  5,  this  produced  a more  gentle  Initial 
rate  of  pressure  rise,  although,  with  the  same  C^,  the  ultimate 
pressure  achieved  was  of  course  the  same. 

As  can  be  seen  from  Fig.  5,  the  combination  of  a parabolic 
dividing  streamline  shape  and  a value  of  C^»0.060  produced  a 
computed  static  pressure  profile  in  good  agreement  with  the 
experiment.  It  should  be  noted,  that  if  a Prandtl  eddy  vis- 
cosity model  is  considered  for  the  shear  layer  region,  a value 
of  C^aO.060  is  equivalent  to  a Prandtl  model  constant  of  0.018, 
which  is  in  good  agreement  with  the  values  reported  for  shear 
layer  analyses. 

Figure  6 shows  a comparison  between  the  computed  and 
measured  centerline  Mach  number  profiles  for  the  same  cold  flow 
case;  it  is  clear  from  these  results  that  somewhat  too  rapid 
mixing  is  predicted  in  the  directed  flow.  Study  of  the  results 
indicates  that  the  problem  may  lie  in  the  modeling  of  the  turbul- 
ent length  scale  distribution.  The  dump  burner  flow  is  real- 
istically a two-scale  problem,  with  one  turbulent  length  scale 
related  to  the  turbulent  shear  layer  region  and  a second  scale 
more  appropriate  for  the  directed  flow.  Thus  for  these  flow- 
fields,  a two-equation  turbulence  model,  in  which  the  length 
scale  distribution  is  computed  through  a transport  equation 
as  part  of  the  solution  may  be  more  appropriate-  and  will  be 
Incorporated  in  future  work. 

Computations  for  a representative  reacting  dump  burner 
flow  have  also  been  carried  out,  again  in  order  to  assess  the 
sensitivity  of  the  results  obtained  to  the  assumptions  used 
in  the  model.  A schematic  of  the  configuration  involved  is  shown 
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in  Figure  7,  which  also  gives  the  conditions  for  the  computa- 
tions. The  stirred  reactor  was  assumed  to  contain  the  products 
of  a fully-reacted,  equilibriiun,  stoichiometric  hydrocarbon- 
air  combustion  process;  the  temperature  of  the  resulting  mixture 
was  arbitrarily  specified.  In  order  to  perform  preliminary 
computations  very  rapidly,  a global  finite-rate  combustion  model 
was  used  in  the  directed  flow.  The  recirculation  zone  geometry 
was  assumed  to  be  the  same  as  was  used  in  the  cold  flow  calculations. 

Figure  8 shows  that  results  of  a representative  computation 
of  a hot  flow  case;  here  a 2200ok  recirculation  zone  temperature 
has  been  assumed.  The  cross-hatching  indicates  the  region  in 
which  significant  combustion  is  occur ing  in  the  directed  flow. 

One  important  observation  made  from  a series  of  such  computations 
such  as  that  shown  here  is  the  importance  of  the  assumed  initial 
shear  layer  thickness  (i.e.,  the  inlet  boundary  layer  thickness) 
and  the  assumed  shear  layer  growth  rate.  The  expression  for 
the  shear  layer  growth  rate  is  Z=a+bx,  where  I is  the  shear 
layer  width  scale;  a and  b are  shown  on  the  figure.  Increasing 
the  initial  thickness  and  growth  rate  of  the  shear  layer  reduces 
the  rate  of  transfer  through  the  shear  layer  and  thus  reduces 
the  temperature  rise  in  the  directed  flow  delaying  the  onset 
of  combustion.  On  the  other  hand,  if  the  assumed  shear  stress 
level  on  the  dividing  streamline  is  increased  (as  in  general 
is  required  for  an  increase  in  shear  layer  growth  rate)  the 
rate  of  momentum  and  energy  transport  is  increased,  and  the 
Increased  residence  time  caused  by  the  resultant  retardation 
of  the  main  flow  allows  more  combustion  to  take  place. 

All  of  the  effects  mentioned  above  manifest  themselves 
quite  strongly  in  the  computed  axial  static  pressure  gradient, 
so  that  this  parameter  can  apparently  be  a tentative  indicator 
of  some  of  the  details  of  the  combustion  process  occur ing  with- 
in the  dump  burner. 
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In  addition  to  the  progress  achieved  in  the  application 
of  the  model  described  above,  during  the  current  contract  period 
close  coordination  has  been  maintained  with  Drs.  F.  D.  Stull 
and  J.  E.  Drewry  at  AFAPL.  This  has  included  providing  a critical 
review  of  the  designs  for  experimental  hardware  currently  under 
test  at  AFAPL,  and  further  recommendations  for  appropriate  in- 
strumentation for  both  hot-and  cold  flow  testing  based  on  the 
results  of  use  of  the  modular  model. 
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VI.  FURTHER  WORK 

Work  on  the  development  and  exploitation  of  the  modular 
combustor  model  will  continue,  with  emphasis  on  the  inclusion 
of  advanced  turbulence  models  to  adequately  assess  the  effects 
of  the  different  turbulence  scales  in  different  parts  of  the 
flow,  as  described  in  the  preceding  section,  and  on  the  inclusion 
of  the  effects  of  swirl  in  the  directed  portion  of  the  flow. 

Exploitation  of  the  modular  model  will  involve  the  use 
of  the  tool  to  develop  a family  of  experiments  relevant  to 
hydrocarbon  fuel  injection  and  combustion,  emphasizing  the 
effects  of  different  fuels  including  methane,  ethylene  and  the 
higher  hydrocarbons,  and  alternate  fuels  derived  from  non- 
petroleiim  sources.  The  focus  of  this  work  will  be  on  the 
generation  of  a data  base  on  steady-state  combustion  performance 
and  flame  stabilization  in  sudden  expansion  burners. 

The  modular  model  will  also  be  used  to  obtain  a preliminary 
assessment  of  the  effects  of  turbulence  levels,  initial  boundary 
layer  thickness,  and  initial  turbulence  scale  distribution  on 
combustor  performance.  The  effects  of  the  degree  of  fuel  vapor- 
ization at  the  combustor  entrance  and  within  the  recirculation 
zones  will  also  be  assessed.  These  results  will  be  used  in 
examining  existing  data  and  identifying  important  new  data 
required  for  model  verification  and  refinement.  In  addition, 
aid  will  be  provided  toward  the  design  of  test  configurations 
and  Instrumentation  (including  probe  design)  based  on  the 
resources  available  at  AFAPL. 

Further,  a number  of  recent  investigations  of  the  use  of 
unified  models  in  combustor  flowflelds  have  pointed  out  the 
need  for  modular  modeling  within  an  overall  elliptic  formulation 
to  overcome  the  grid  resolution  problems  commonly  encountered 
in  the  development  of  unified  codes  for  reacting  flowflelds. 

Thus  the  development  of  a unified  model  specifically  designed 
for  the  sudden-expansion  combustor  geometry  which  incorporates 
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modular  features  represents  a natural  development  and  extension 
of  our  current  model  formulation. 

Finally,  work  done  to  date  indicates  that  additional  detail 
is  required  in  modeling  the  fuel  injection  and  vaporization  process. 
Thus,  modular  model  computations  will  be  utilized  to  develop  a 
more  detailed  characterization  of  the  fuel  spray  process. 
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